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We investigate the angular distribution of the lepton pair in the process pp ^ "f* + 
X +X, where the virtual photon is produced at high transverse momentum. 

The angular distribution of the leptons is very sensitive to possible nonperturbative 
effects, such as a nontrivial vacuum structure of QCD, and offers a good chance 
to test such effects. We present complete ©(a^) calculations of the decay lepton 
distributions in the lepton pair rest frame. An order 0{as) Monte Carlo study of the 
lepton angular distributions, with acceptance cuts and energy resolution smearing 
applied to the leptons, is also presented. 
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I. INTRODUCTION 



The production of a virtual photon at hadron coUiders, along with the subsequent decay 
into a lepton pair, provides a unique opportunity for testing perturbative Quantum Chro- 
modynamics (QCD). In particular, the measurement of the angular distribution of the decay 
leptons provides a detailed test of the production mechanism of the virtual photon. Exper- 
imental studies of the angular distribution of the lepton pair have been made by the NAlO 
collaboration at CERN [Q] and the Chicago-Iowa-Princeton collaboration at Fermilab 0. 
The data from these experiments do not agree with the predictions of the QCD improved 
parton model. An investigation of the decay lepton distribution at the Tevatron center of 
mass energy would certainly be interesting in view of these results. In this paper, we present 
next-to- leading order QCD predictions for the angular distributions of the lepton pair in the 
process pp —>■'-/*+ X l^l" -|- X at Tevatron energies together with a leading order Monte 
Carlo study in which acceptance cuts and energy resolution smearing are applied to the 
leptons. 

A basic assumption in the calculation of the Drell-Yan process is the factorization 
hypothesis, which asserts that the hadronic cross section can be calculated by convoluting the 
cross section for the corresponding parton level reaction, gg 7* ^ with the parton 

distribution functions of the initial state hadrons. Great theoretical efforts have been made to 
prove this factorization assumption within perturbative QCD [Q. However, in Ref. it has 
been argued that nonperturbative QCD effects may spoil the factorization. Nonperturbative 
vacuum fluctuations could generate a background field (like the color domains described in 
the first paper of Ref. f^) which could induce spin and color correlations between the partons 
of the incoming hadrons. These effects need not vanish in the high energy limit, i.e., they 
can be "leading twist" effects (for reviews see Ref. 0). 

It has recently been shown that the decay lepton distribution of a vector boson 
produced in hadronic collisions is very sensitive to nonperturbative effects and offers a good 
chance to test the factorization hypothesis. As mentioned earlier, the experimental results 
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in Refs. [Q and for the decay lepton distribution from a virtual photon produced in pion- 
nucleon scattering are not in agreement with next-to-leading order QCD predictions [0 and 
therefore give an indication that factorization may be violated. In Ref. [0 it has been shown 
that these results can be explained by transverse spin correlations between the initial state 
partons, which may be induced by the nontrivial vacuum structure in QCD. Pion bound state 
effects have been discussed in Ref. [§] as another possible explanation for the discrepancy. 

It would clearly be interesting to check whether the decay lepton distributions from gauge 
bosons produced at the Tevatron are in agreement with the factorization assumption. To 
do this test, reliable theoretical predictions for the decay distributions are mandatory. In 
this paper, we present complete next-to-leading order predictions for the angular 

distributions of the leptons produced via the decay of a high virtual photon, i. e., the Drell- 
Yan process ^ 7* + X — + X is calculated to C(a^). The calculation is based on 
the assumption of factorization, i.e., the calculation is done in the standard QCD-improved 
parton model. When the virtual photon is produced with no transverse momentum, the 
zero order Drell-Yan subprocess, gg 7* ^ predicts a 1 + cos^ 6 distribution for 

the leptons, where 6 is the scattering angle in the parton center of mass frame (the 7* rest 
frame). For virtual photons produced with transverse momentum (balanced by additional 
gluons or quarks), the event plane spanned by the beam and virtual photon momentum 
directions provides a convenient reference plane for studying the angular distributions of the 
decay leptons. The angular distribution now has the general form 

, ^ ~ (l + cos02) + -Ao(l-3cos2^) + Ai sin 2^ cos + - sin20cos20, (1) 

d(p d cos 9 2 2 

where 6 and denote the polar and azimuthal angles of the decay leptons in the virtual pho- 
ton rest frame. The coefficients Ai are functions of the transverse momentum and rapidity of 
the virtual photon (measured in the laboratory frame) and vanish in the limit Pt{i*) ~^ 0. 

The remainder of this paper is organized as follows. In Sec. II we discuss the formalism 
for describing the angular distributions of the decay leptons, give an overview of the 0{al) 
calculation of the process — 7* + X — ^ i^i" + X, and define two choices for the z- 
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axis in the lepton pair rest frame. In Sec. Ill we present 0{al) numerical results for the 
coefficients Ai and show that the QCD corrections to the coefficients are fairly small. This 
is because the coefficients Ai are ratios of helicity cross sections [see Eq. (|Tl])] and the QCD 
corrections tend to cancel in these ratios. A leading order Monte Carlo study of the 

decay lepton distributions for the Tevatron center of mass energy, with typical acceptance 
cuts and energy resolution smearing applied to the leptons, is also given. Conclusions and 
summary are given in Sec. IV. Finally, there is an appendix containing next-to-leading order 
matrix elements which contribute to virtual photon production but have not appeared in 
the literature. 

II. FORMALISM 

The formalism used in our calculations is discussed in this section. First the methodology 
for describing the angular distributions of the leptons is reviewed, then an overview of the 
0{al) calculation of the process — 7* + X — s> + X is given, and finally, two choices 
for the z-axis in the lepton pair rest frame are defined. 

A. Angular distributions 

We consider the angular distributions of the leptons coming from the leptonic decay 
of a virtual photon produced with non-zero transverse momentum in high energy proton- 
antiproton collisions. For definiteness we take 

p(Pi) + p{P2) 7*(Q) + ^ ^ ^+(^1) + r (4) +X, (2) 

where the quantities in parentheses denote the four-momenta of the particles. At leading 
order the parton subprocesses 

q + q^'y* + g , q + g^i* + q, (3) 

contribute to high p-^ virtual photon production. 
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In the parton model the hadronic cross section is obtained by folding the hard scattering 
parton level cross section with the respective parton densities: 



dQ^tllyd^* =^,h''' dQ^t2dn* (^^^^' "^^^^^ 



(4) 



where the sum is over a,b = q, q, g. fl^{x, fi'j^) is the probability density for finding parton a 
with momentum fraction x in hadron h when it is probed at scale fi'j^. The parton level cross 
section for the subprocesses in Eq. (|) are denoted by ddat- Note that in the framework 
of the parton model the incoming partons are assumed to be unpolarized in spin and color 
(for unpolarized initial state hadrons). Furthermore, one neglects the transverse momenta 
of the incoming partons. (This is in contrast to the model discussed in Ref. where it is 
assumed that a nontrivial structure of the QCD vacuum correlates the spins and momenta 
of the incoming partons. Given such a correlation, it has been shown that these correlations 
can drastically affect the polarization of the produced boson ||^.) 

Denoting hadron level and parton level quantities by upper and lower case characters, 
respectively, the hadron and parton level Mandelstam variables are defined by 



and 



S=iPi + P2f , T = 


{Pi - Qf 


U={P2- Qf , 


(5) 


S = (pi +P2f 


= X1X2S , 






t = {Pi- Q? 


= Xi(T- 


Q^) + Q^ , 


(6) 


U = {P2- Q? 


= X2{U- 


Q') + Q' , 





where pi = XiPi and p2 = X2P2- The rapidity y of the virtual photon in the laboratory 
frame can be written 

and the transverse momentum pj, of the virtual photon is related to the Mandelstam variables 
via 



The angles 9 and (j) in dVl* = d cos 6 d(f) are the polar and azimuthal decay angles of the 
leptons in the virtual photon rest frame, measured with respect to a coordinate system to 
be described later. The angular dependance of the differential cross section in Eq. (§) can 
be written (see Ref. Q for details) 

IStt da _ da^^'^ ^ 6') + ^'^^ (13 ^ ^) (9) 

3 dQ"^ dp^ dy d cos dcj) dQ^ dp^ dy dQ"^ dp\ dy 

+ doZv^ dy ^ ^ ^'^ + rfO W dy • 
The unpolarized differential production cross section is denoted by a'^'^^ whereas o"^'^'^ 
characterize the polarization of the virtual photon, e.g., the cross section for longitudinally 
polarized virtual photons is denoted by a^, the transverse interference cross section by 
0"^, and the transverse-longitudinal interference cross section by (all with respect to the 
chosen z-axis of the lepton pair rest frame). The hadronic helicity cross sections ^Q2^pi 
in Eq. (P) are obtained by convoluting the partonic helicity cross sections with the parton 
densities: 

rfo"" r , , , , , „v.„, sdcr 



Introducing the standard angular coefficients |11 



. _2da^ _2V2 da^ . _ ^ da^ . . 

^0-^^, ^i-^^cTT^' ^2-^^, (11) 

the angular distribution in Eq. (RI) can be conveniently written 



da 3 da^+^ 



dQ"^ dp^ dy d cos 6 dcj) IGvr dQ"^ dp^ dy 



:i + cos2 0)+ ^Ao (l-3cos2^) (12) 



+ Ai sin 2^^ cos + -A2 sin^^^cos20 



Integrating the angular distribution in Eq. (|T2D over the azimuthal angle yields 



where 



3 rf(T^+^ 

O — — 



Integrating Eq. 



8 rfQ2 ^p'^ 

over the polar angle 6 yields 
da 1 da^+^ 



a 



2-3Ao 
2 + Ao 



dQ"^ dp\ dy d(j) 27r dQ'^ dpj^ dy 



[1 + (3cos2(f)) 



where 



A, 



(14) 



(15) 



(16) 



By taking moments with respect to an appropriate product of trigonometric functions it is 



possible to disentangle the coefficients Ai. A convenient definition of the moments is |[T3 

/ da{prp,y,9,(j)) m d cos 9 dcj) 



J da{prp,y,9,(f)) dcos9i 
which leads to the following results: 



(1) = 1 



(-(l-3cos^^^)) 



(sin 29 COS0) 
(sin^ 9 cos 20) 



— Mo — 

20 V 3 



(17) 

(18) 
(19) 
(20) 
(21) 



B. Next-to-leading order cross section 



At C(as) the following partonic tree level and one-loop processes contribute to the par- 



tonic helicity cross sections 



dQ^ dt du 



in Eq. (|T0|). 



tree level contributions : q + q 

q + q 

Q + g 

q + q 
9 + 9 



l* + 9 + 9, 
l* + q + q, 

1* + q + 9 , 

1* + q + q, 
1* + q + q, 



(22) 



one-loop contributions : q + a 7* + g , 

(23) 

q + g 1* + q. 

The second-order contributions in Eq. ( pSf ) come from the interference of the one-loop dia- 
grams with the leading-order diagrams. Let us briefly sketch the technical ingredients that 
go into our calculation (for more details, see Ref. |^). For the C(a^) tree level contributions 
in Eq. (|2^) we introduce the variable 

S2 = ih + kif = {pi+p2-Qf = s + t + u-Q\ (24) 

in addition to s, t, u defined in Eq. (P); S2 is the invariant mass of the system recoiling against 
the virtual photon. 

To obtain the pj. distribution of the virtual photon, the 0[a1) tree level diagrams must 
to be integrated over the phase space of the two final state partons with the pj. of the virtual 
photon held fixed. The integration over the recoiling partons is most easily performed in the 
(^1^2) center of mass system by integrating over the solid angle dVL^^k^. Since all partons are 
massless, coUinear divergencies appear after integrating over dflkik2- Soft gluon singularities 
show up as poles in the variable S2- A finite next-to- leading order (NLO) partonic helicity 
cross section is derived in the following manner. 

• Infrared and collinear divergencies associated with final state partons cancel among 
loop and tree diagrams. 

• Collinear initial state divergencies are absorbed into the parton densities, i.e., they are 
removed by renormalizing the parton densities, which introduces a factorization scale 
dependence into the parton densities /(x,/i|,). 

• Ultraviolet divergencies are removed by MS (Modified Minimal Subtraction |p!0[ ) 
renormalization, which introduces a renormalization scale dependence into the strong 
coupling constant ^^(/i^). 

8 



Following Ref. 0, we introduce the following notation to list the partonic helicity cross 
sections ^^jg^ in Eq. (|10D [a e {U + L, L, T, I}]. 
0{as) Born contributions: 

, a.Born ry-y* 

_^.^„.^(. + < + „_Q^)r^(B). (25) 
virtual corrections: 

iQ^u = ^ ^ + * + « - «'^) Vi^) . (26) 

0{al) tree graph corrections: 

The subscript ab stands for the initial state parton pair, i.e., a and b denote a quark, anti- 
quark, or gluon. The initial state coUinear singularities have been factorized from Eq. (pTj) 
at a scale /i^. The constants Kjf^ and V{e) are given by (we work in n = 4 — 2^ dimensions) 



. ^ Ifor Cf (1-£) /fir^y /fO^V ,„„, 

«« 3 87rQ2 Afo r(l -e) I M tii ' ' ^ ' 



T^(B,V,T) are the partonic helicity matrix elements for the Born, virtual, and tree contribu- 
tions. The explicit form of the partonic helicity matrix elements depends on the choice of the 
z-axis in the 7* rest frame. Covariant projection cross sections da'^ [P G {U+L, Li, L2, L12}] 
for W boson production have been calculated to in in Ref. 0. From these results one 

can obtain the relevant helicity cross sections in Eqs. (^5D-(^) for any given 7* rest frame 
[see Eqs. ( |Al| ) and (|A3|) ]. However, there are some further contributions in the case of 7* 
production, which do not contribute to W boson production because of charge conservation. 



The diagrams can be found in Figs. 8 and 9 of Ref. ||^. The remaining analytical projec- 
tions for these contributions (after integration over dQf^-^k^) are listed here in an appendix. 
The folding of the NLO parton level cross sections with the respective parton densities is 
straightforward and we do not list the combinations here. 



C. The z-axis in the lepton pair rest frame 

Before we present numerical results for the angular distributions, it is necessary to discuss 
the choice of the z-axis in the lepton-pair rest frame. We will discuss two different choices: 
the CoUins-Soper (SC) frame and the Gottfried- Jackson (GJ) frame [Q. In the CS 
frame the z-axis bisects the angle between Pi and — P2: 



Pi = El (sin 7^5, 0, cos Jcs) ) 

— * 

P2 = E2 (sin 7(^5, 0, - cos 7^5) , 



with 



with 



(32) 



sin 7cs = - V ^ ~ ^cs > (34) 

In the GJ frame the z-axis is chosen parallel to the beam axis: 

Pi = El (0, 0, 1) , 
P2 = E2 (sin7Gj, 0, cos7gj) , 



(36) 



^^-^ " {T-Q^){U-Q^)- pI + Q^' ^^^^ 



sin 7g J = - V ^ ~ J , (38) 
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Note that the CS and G J frames are related by a rotation about the y-axis. In the laboratory 
frame, the ^-direction is defined by the proton momentum and the x-direction is defined by 
the transverse momentum of the virtual photon. 



III. NUMERICAL RESULTS 

In this section numerical results are presented for high production and leptonic decay 
of a virtual photon at the Tevatron collider center of mass energy (v^ = 1.8 TeV). The 



numerical results have been obtained using the MRS set A |Tj] parton distribution functions 
with = 230 MeV. For our NLO predictions, we use the two-loop formula for with 
five favors. If not stated otherwise, the renormalization scale fij^ and the factorization scale 
in Eq. (P have been taken to be /ijj = /i^ = [Q^ +Pt(7*)]/27 where Q and ^^(7*) are 
the invariant mass and transverse momentum, respectively, of the virtual photon. We work 



in the MS factorization scheme. 

We begin with numerical results for the coefficients Ai, a, and (3 in Eqs. (|12|), (p!3[), 
and (p!5D. Figure la shows the coefficients Ao,Ai, and A2 in the CS frame as functions of 
Pxi'-f*) for an invariant mass Q of the photon fixed to Q = 10 GeV. The reason for choosing 
this fairly large invariant mass is to minimize the effect of the acceptance cuts which will 
later be imposed on the decay leptons [see below]. The dotted lines are leading order (LO) 
predictions and the solid lines are next-to-leading order [(^(ag)] predictions. The 
coefficients Aq and A2 are increasing functions of Pt{i*) and the deviations from the zero 
order expectation [Aq = = at ^^(7*) = 0] are quite large, even at modest 



values of pj,['~)*). It has been noted in Ref. |T5| that the coefficients Aq and A2 are exactly 



equal at LO (dotted line). This is no longer true at NLO, but the corrections are fairly 
small, especially the corrections to Aq. The 0{a1) corrections to A2 are negative and about 
20% the size of the LO result. Note that Aq originates from the longitudinal polarization of 
the virtual photon, whereas A2 receives contributions from a transversely polarized virtual 
photon [all with respect to the z-axis of the chosen lepton pair rest frame]. The deviation 
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of Ai (interference of longitudinal and transverse virtual photon polarizations) from zero is 
small in the CS frame, even at large values of ^^(7*). 

Figure lb shows numerical results for the coefficients a and /5 [see Eqs. (|T3|) and (|1^)] as 
a function of ^2^(7*). The coefficients are again very sensitive to the transverse momentum 
of the virtual photon and the deviations from the zero order predictions (a = 1 and /3 = 0) 
are again large. The NLO corrections to a are small over the whole range of ^^(7*). 

Numerical results for the coefficients Ai, a, and f3 in the GJ frame are shown in Fig. 2. 
The corrections to Aq and a are larger than in the CS frame, however, in both frames the 
corrections do not dramatically change the LO results. This is because the coefficients Ai are 
ratios of helicity cross sections [see Eq. ([TT|)] and the large QCD corrections in the individual 
helicity cross sections tend to cancel in the ratios. This can be seen from Fig. 3, where the 
i^-factors for the helicity cross sections and are shown as functions of ^^(7*) 

in the CS frame (Fig. 3a) and GJ frame (Fig. 3b). The i^-factor is defined as the ratio of the 
NLO [0(a^)] differential cross section to the LO [0(a^)] differential cross section. Results 
are shown in Fig. 3 for two different choices of the scale = = /i^: the upper curves 
correspond to /i^ = 1/4 [pt{i*) + Q"^] and the lower curves correspond to /i^ = 4 [pt{i*) + Q'^]- 
The invariant mass of the virtual photon has again been fixed to Q = 10 GeV. The ii"-factor 
for a^'^^ and ranges from 0.9 to 1.6 in the CS frame depending on ^^(7*) and the choice 
of the renormalization and factorization scale. However, the i^'-factor is almost the same 
for a^~^^ and a^, thus Aq, which is proportional to the ratio /a^^^, is not effected by 
the corrections. In the GJ frame, the i^-factors for and cr^^^ differ more, and thus 
the corrections to Aq are larger in the GJ frame (see Fig. 2b). The K-factors for o"^ are 
particularly different from the ii'-factors for a^^^ , which explains the large deviation of A2 
from the LO result Aq = A2. Note that the large i^'-factors are due to large logarithms [like 
ln(s/(5^)] in the NLO matrix elements The i^-factors decrease with increasing invariant 
mass of the virtual photon. 

To give a feeling for the numerical contributions from the different partonic subprocesses. 
Fig. 4 shows the fractional contributions of the partonic subprocesses to and a 
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for 7* production as a function of Pt(7*) for Q = 10 GeV. The curves labeled A and B 
correspond to the LO results from the subprocesses qq '~f*g and qg '~f*q, respectively. 
The curves labeled C through G correspond to the 0{al) contributions from the subprocesses 
listed in Eqs. ( p2D and p^). Note that the NLO contribution from the qg initiated subprocess 
(curve E) is more important than the LO contribution from the qq subprocess (curve A) at 
large values of pj,(7*). The other O(a^) subprocesses are fairly small and even give negative 
contributions at small values of ^7^(7*). This is due to the factorization of the coUinear initial 
state singularities which dominate these subprocesses at low values of Pxil*)- Although the 
total contributions from the 0{al) subprocesses are large, they do not have a dramatic 
effect on the decay lepton distribution (see Figs. 1 and 2), i.e., the polarization of the 
virtual photon as a function of its transverse momentum is almost unchanged by higher 
order corrections. 

We want to point out that direct measurements of the coefficients Ai for virtual photon 
production in pion-nucleon scattering at \^ ~ 19 and 23 GeV are not in agreement with 
the LO QCD predictions The NLO corrections to the coefficients are very small in 

this case (see Fig. 7 in Ref. [0) and the results in Refs. [|^ and 0] cannot be explained by 
standard Drell-Yan production supplemented with QCD corrections. 

We now turn our attention to the cos^ and (p distributions of the decay leptons. Since 
the effects of the NLO corrections are small (in particular for the cos 6* distribution in the 
CS frame; see the coefficients Aq and a in Fig. 1), it is sufficient to use LO matrix elements 
in our Monte Carlo study of the lepton decay distributions. Results will be shown for the 
and cos 6 distributions of leptons originating from the decay of a virtual photon produced 
with finite transverse momentum in pp collisions at the Tevatron center of mass energy. To 
demonstrate the effects of cuts, results are shown first without cuts and then with typical 
acceptance cuts imposed on the leptons. These cuts are necessary due to the finite acceptance 
of the detector. 

Measurement uncertainties, due to the finite energy resolution of the detector, have been 
simulated in our calculation by Gaussian smearing of the lepton four-momentum vectors 
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with standard deviation cr. The numerical results presented here were made using a values 
based on the CDF specifications ||16| . The energy resolution smearing has a negligible efi^ect 
on the (f) and cos 9 distributions. 

Figure 5 shows the normalized (f) and cos^ distributions of the decay leptons from 7* 
production for three bins in the transverse momentum of the virtual photon. The invariant 
mass of the photon has been integrated over the range 10 GeV < Q < 12 GeV. No cuts 
or smearing have been applied to the results in the figure. The curves in Fig. 5 can be 
obtained from a and j3 in Fig. lb. [Note, however, that Fig. 1 shows the coefficients a 
and (3 for a fixed value of Q = 10 GeV.] For example, since (3 is an increasing function of 
Pt(7*); amplitude of the distribution increases withpj,(7*). Likewise, a is a decreasing 
function of ^^.(7*), starting out positive and ending up negative, thus the curvature of the 
cos^ distribution in Fig. 5b is positive for the two lowest Pxil*) bins and is nearly zero for 
the Pt(7*) > 6 GeV bin. Note that if the virtual photon was to decay isotropically, the 
and cos 6' distributions would both be flat. 

The effect of acceptance cuts on the angular distributions is illustrated in Fig. 6 which 
shows the and cos6 distributions of the decay leptons for the same bins in ^^(7*) as in 
Fig. 5, but now with energy resolution smearing and the cuts 

> 2 GeV, \y{e)\ < 2.5 , \y{j*)\ < 1-0 • (40) 

The cuts on the leptons are necessary due to the finite acceptance of the detector. The photon 
rapidity cut has been imposed because we find that polarization effects are highlighted when 
the virtual photon is in the central rapidity region. The cuts on the leptons, in particular 
the Pt(^) cut, have a dramatic effect on the shapes of the distributions. The shapes of the 
angular distributions are now governed by the kinematics of the surviving events. Only 20% 
of the events pass these cuts. The cuts in Eq. (^0]), which are applied in the laboratory 
frame, introduce a strong (p dependence. The "kinematical" distribution in Fig. 6a is 
very different from the "dynamical" distribution in Fig. 5a, in fact, the peaks and valleys 
are interchanged in the two distributions. The only remaining vestiges of the polarization 
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effects in the distribution are the dips in the high Pt{i*) curve (sohd hne) at = 90° and 
270°. The cos 6* distributions with cuts in Fig. 6b are also different from the corresponding 
results without cuts in Fig. 5b, in particular for |cos^^| ^ 0.5. However, for |cos0| < 0.5 
and small Pt{i*) (dotted and dashed curves in Fig. 6b) the polarization effects in the cos^ 
distributions are still visible. As stated earlier, the effect of energy resolution smearing is 
negligible; the drastic changes in the shapes of the distributions are due to the cuts, especially 
the Pt(^) cut. The polarization effects diminish as the invariant mass of the virtual photon 
decreases. The LO cross section (summed over i = e,fi) for the three Pxil*) bins, 2 GeV 
< Pt(7*) < 4 GeV, 4 GeV < ^^(7*) < 6 GeV, ^^(7*) > 6 GeV, are 50 pb, 20 pb, 26 pb, 
respectively. However, these numbers should be multiplied by the i^-factor for cr'^+^ as 
shown in Fig. 3. Furthermore, soft gluon resummation effects will also be important for the 
production cross section at low ^^.(7*). 

The cuts in Eq. (|40| ) are actually quite weak. Figure 7 shows the and cos 6 distributions 
with the stronger and more realistic cuts 

p^(£)>5GeV, W)\<1.0, |2/(7*)|<1.0, (41) 

for the same three ^^(7*) bins as in Fig. 5. The curves are now very different from the 
curves in Fig. 5 and no traces of polarization effects are left in Fig. 7. 

We have also analyzed the effect of the cuts by using the correct matrix element for 7* 
production, but with isotropic decay of the virtual photon, i.e., neglecting spin correlations 
between 7* production and decay. The angular distributions in this case are similar to the 
ones shown in Figs. 6 and 7 for the full matrix element; the remnant polarization effects 
discussed in Fig. 6 are of course absent. 

In Fig. 8 we show ratios of the and cos6 distributions for the same bins in Pj'{'j*) 
as in Fig. 5; the distribution with full polarization has been divided by the corresponding 
distribution obtained with isotropic decay of the virtual photon. Cuts and smearing are 
included in both cases. The large effects from the cuts are expected to almost cancel in this 
ratio. In fact, we nearly recover the and cos 6 dependence of Fig. 5 which contains no cuts. 
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Fig. 9 shows the ratio of the and cos 6' distributions with full polarization to the 
corresponding distributions obtained with isotropic leptonic decay for the virtual photon 
for the cuts in Eq. (pj) for the high Pxil*) bin, i.e., Pxil*) > 6 GeV. The kinematical 
effects in the two low PTil*) bins are very large and the cuts remove all events around 
= 0°, 180°, 360° as well as for large cos 6' values [see Fig. 7], thus it is impractical to form 
ratios for the two low p^il*) bins. However, the ratios for prp{'~f*) > 6 GeV [see Fig. 9] 
once again contain most of the polarization dependence seen in the solid curves of Fig. 5. 
The additional dips in the distribution at = 0°, 180°, and 360° in Fig. 9a are due to the 
kinematical cuts. Thus even in the presence of large acceptance cuts it may be possible to 
highlight the polarization effects in the experimental results by dividing the experimental 
distributions by the corresponding Monte Carlo distributions obtained using isotropic 7* 
decay. 



IV. SUMMARY 

The polar and azimuthal angular distributions of the lepton pair arising from the decay 
of a virtual photon produced at high transverse momentum in hadronic collisions have 
been discussed. In the absence of cuts on the final state leptons, the general structure 
of the lepton angular distribution in the virtual photon rest frame is determined by the 
polarization of the virtual photon. In perturbative QCD, the structure is described by four 
helicity cross sections, which are functions of the transverse momentum and rapidity of 
the virtual photon. We have calculated to Ola^) the angular coefficients which govern the 
lepton angular distributions and find that the corrections are relatively small in both the 
CS and GJ frames, especially in the CS frame. This is because the angular coefficients Ai 
are ratios of helicity cross sections and the large QCD corrections in the individual helicity 
cross sections tend to cancel in the ratios. 

We have also studied the angular distributions of the leptonic decay products of a high 
Pj. virtual photon when acceptance cuts and energy resolution smearing are applied to the 
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leptons. When acceptance cuts are imposed on the leptons, the shapes of the lepton angular 
distributions are dominated by kinematic effects and the residual dynamical effects from 
the virtual photon polarization are small. The kinematic effects become more dominate as 
the cuts become more stringent and as the invariant mass of the photon decreases. Energy 
resolution smearing has a negligible effect on the angular distributions. 

Polarization effects can be maximized by minimizing the cuts on the decay leptons, 
however, this strategy is severly limited since cuts are needed due to the finite acceptance 
of a detector. Polarization effects are also more pronounced when the virtual photon is in 
the central rapidity region. Alternatively, it may be possible to highlight virtual photon 
polarization effects by "dividing out" the kinematic effects, i.e., if the histogrammed data 
is divided by the theoretical result for isotropic virtual photon decay, the resulting ratio is 
more sensitive to polarization effects. 
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FIGURE CAPTIONS 

Fig. 1 Angular coefficients for 7* production and decay in the CS frame as a function of the 
transverse momentum PxiY) fo^" Q = 10 GeV and =1.8 TeV. Part a) shows the 
angular coefficients Ao,Ai, and A2 and part b) shows the angular coefficients a and 
The dotted lines are LO predictions {Aq = Ai at LO) and the solid lines are NLO 
predictions. No cuts or smearing have been applied. 

Fig. 2 Same as Fig. 1 but for the GJ frame. 

Fig. 3 i^-factors for the helicity cross sections da'^~^^, da^, and da^ for 7* production and 
decay as a function of the transverse momentum Pr(7*) for Q — 10 GeV and — 
1.8 TeV. Parts a) and b) are for the CS and GJ frames, respectively. Results are 
shown for two different choices of the scale = fi^p = the upper lines are for 
/x^ = 1/4 \pt{j*) + Q'^] and the lower lines are for //^ = 4 [pt{i*) + Q^]- 

Fig. 4 a) Fractional contributions to da^^^ for 7* production in the CS frame as a function 
of the transverse momentum ^^(7*) for Q — 10 GeV, \fS — 1.8 TeV, and jj? — 
1/2 \pt{i*) + The curves are labeled according to the contributing subprocesses. 

The LO subprocesses are: 
(A) qq Yg and (B) qg YQ- 

The NLO subprocesses are: 

(C+D) {qq Yg) + {qq ^ Ygg) + {qq ^ Yqq), 

(E) {qg ^l*q) + {qg ^l*qg), 

(F) qq YQQ, 

(G) gg YqQ- 

b) Same as a) but for da^. 

c) Same as a) but for da^. 

Fig. 5 a) Normalized and b) normalized cos 9 distributions of the leptons from 7* decay in 
the CS frame with 10 GeV < Q < 12 GeV. Results are shown for three bins in ^^(7*): 



2 GeV < Pt{Y) < 4 GeV (dotted), 
4 GeV < pr(7*) < 6 GeV (dashed), 
Pt{Y) > 6 GeV (solid). 

No cuts or smearing have been applied to the decay leptons. 
Fig. 6 Same as Fig. 5 but with smearing and the cuts p-^ii) > 2 GeV, \y{i)\ < 2.5, and 

|y(7*)l < 1-0. 

Fig. 7 Same as Fig. 5 but with smearing and the cuts p^ii) > 5 GeV, \y{i)\ < 1.0, and 
\y{Y)\<1.0. 

Fig. 8 Ratios of distributions in the CS frame obtained with full polarization to those obtained 
with isotropic decay of the 7* for the same Pxil*) bins as in Fig. 5. Parts a) and b) are 
the ratios for the and cos 9 distributions, respectively. Energy resolution smearing 
and the cuts p^ii) > 2 GeV, \y{e)\ < 2.5, and \y{Y)\ < 1-0 are included. 

Fig. 9 Same as Fig. 8 but with the cuts p^ii) > 5 GeV, \y{e)\ < 1.0, and |y(7*)| < 1.0. Only 
the bin for Pr(7*) > 6 GeV is shown. 
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APPENDIX A: NLO MATRIX ELEMENTS 



In this appendix, we present the remaining NLO matrix elements which do not contribute 
to W boson production and are thus not hsted in Ref. [Q. It is convenient to calculate 
covariant projection cross sections da^ [f3 G {U + L, Li, L2, L12}] from which one can deduce 
the helicity cross sections da°' [a G {U + L, L, T, I}] in Eqs. (p5D-(^7D for any given 7* rest 
frame by a transformation matrix {M)ai3- For the CS frame one has 



I 



da^ 
da^ 



\ 



I 











1 







4 003"^ '^CS 

1 (1+cos^ ics) 

2 8 sin"* 'yes cos^ -y^g 







-1 



4 COS'' 7CS 4 cos^ "/CS 

(1+cos^ JCS) (l-3cos^7cg) 

— — n T: : — ^ o 



CS 



sin 7CS cos fcs 4\/2 sin 7cs cos 7cs 



8 sin"^ yes cos^ 7^5 8 sin'' 7c7g cos^ ycs 





-1 



(Al) 



where 



cos 7cs 



sm 7CS- 



The results for the GJ frame can be obtained from: 



1 — COS^ 7c5 . 



(A2) 



da^ 
da^ 
da^ 



\ 




1 













1 (l+cos^7Gj) 

2 2sin^7Gj 

- cos 7G J 



1 



COS7GJ 
sin^ 7GJ sin^ yaj 

n 1 



GJ 



V2 



sm7Gj 



2\/2 



sm 7GJ 



V 



(A3) 



and 



cos 7g J 



sm 7g J 



1 - cos2 7g J . 



(A4) 



(t-g2)(«-g2) ' 

All of the partonic projection matrix elements T^^ are listed in Ref. ||^, however, there are 
some interference contributions which do not contribute to W boson production because of 
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charge conservation. The analytical results for these additional partonic projection matrix 
elements are listed here. 



I) Diagrams 2(^5 + F6)*(F7 + Fg) for qq YQQ 

The diagrams are shown in Fig. 8 of Ref . . The result for the interference of these diagrams 
differs for vector-vector (relevant for 7* production) and axial-axial (relevant for Z boson 
production) couplings. They do not contribute to W boson production. The vector-vector 
coupling contribution for the projection matrix elements Tgg for these diagrams are denoted 
by D^^yy. All quantities in the following formulae are defined in appendix F of Ref. 0. 



Dl!,\\,{s,t,u,Q'') 



S(So — Ut) U + t 



ut 



dfdu - ^^r— ^ - {s + S2 - 



2ut 



2utX^ 



h 
A 



(3s + 2s2 



u + t 



ftuQ'^ 

ut 



2ut 

[2{S2 - t) 



2{s-Q^)ds - s{u + t) 



[n-tf 
2utX^ 



' + sl + s' 



— fxt 

+ fst 



t 

dstds 
ut 



S2 - Q) + (s2 - ty + si + s' 



(2Q2 _ t) (e + 2{sl + s^- S2t) 



u 



u 



D^jyyis, t, u, Q') = (^) ( ^ { ^d,td, [-48s^ + 8s^{13s2 -7t- 4u) 

- As^ (3s2(8s2 - 9t) - 13u(s2 - t) - 2u^ + St^) 

- 2s (2u(20s^ - 27s2t + 8t^) - (18s^ - 21s2t + 4t^)(2s2 - t) - (9s2 - 8t) 
- (5m - 8s2 + 3t)(s2 - m)(s2 -t){u- 4S2 + 3t) 

+ iff '"^ [4s^(m + 3t - 2S2) + 2s (7s2(m - t) + 3(^2 - u^) - 2{s\ - ut) 

+ (2S2 -U-t) (7(S2 - Uf - (S2 - tf) 

+ ij;'''^ (S2 - u) [2 {s[u - 2s2) + 2(2s2 -u-t){s2- u)) 
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+ dstu (4s^ - 2s{u + i - S2) + (s2 -u){u + S2- 2t)^ 
+ Hi^'^^ (s2 - uf [(2s2 -u-t)+ udsdst (-4s^ - 2s{3u - 4s2 + t) 

+ (S2 - t){3u - 4S2 + t))] + i/f '^^ M4t(s2 - Uf 



+ i/^^'^^ 2(s2 - t) [(4s2 - 2m - 3t)s - {S2 - t){2s2 -u-t) 
+ (l6s2 + 2s(3ii - lls2 + 7i) + 2u{u - 7s2) + t{10u + 1 - I2S2) + l^sj) 
+ H^^'^^ {s2 - tf [-{2s2 -u-t)+ udstds (-I2s^ - 2s{bu - 12s2 + 7t) 

- {3U - AS2 + t){u- 3S2 + 2t))] + H^2'^'^ u{s2 - tfdst 
+ 2dt [4.8^82 + 2s (s2(2k - 3i - S2) + 2t^) 

- {S2 - t) {{t - S2f - 3{u - S2f)\ + dst^ (24s2 - 4s(4s2 - 3t) 

+ {U - tf + IQ{S2 - U){S2 - t) 

- ^ I ftu^dsudst {{S2 - U){S2 -t) - SS2) {S -U + S2) 
+ fsuUs{s -U + S2)dsdst - fxudsu{s - U + S2){S2 - u + dgUs) 

- fstdsuS {u-s -S2 + dsu{s2 - u - 3s)) 
+ fxtdsdstU ((g' - S2){S2 -t)-s^)^\ 



D^iyyis, t, U, Q') = D^iyyis, U, t, Q') 



D^fwis, t, u, Q') = (i) ( [ - ds [%s\Zu' - hut - 6US2) 

+ 4s^ (iu^iSu - 2t) - l?>us2{?>u -t) + ZQusfj 

+ 2s^ {2u^{Su + 9i) - u'^{l?,us2 + 22^^) + 5^282(9^ + 2OS2) - l2usl{At + 3s2)) 
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+ s (3it^(3w + t)- Au^{13us2 + 19t^) 

+ 12uHs2{l7t + 8u) + 8x^^52(11^* - 52i) - Ausl{7u - 49i + 6S2)) 

+ u\u + 2t- 7s2) - 3u'^{9t^ - 4:sl) + 5uHs2{5u + 46t) + 10u^sl{u - 4s2) 

+ 8us^{3s2 - Ut) + 2uHsl{179s2 - 129t - 93m) - 32uH^ 

+ h[^'^^ 2{s2 - u) [ 3s'^u{u -3t- 2S2) + s (Au^u - 3^2) 

- 3ut{u + t - 4s2) + 4:sl{2u -t) + 2^82) + u^{u + t - AS2) 

+ 2ut{2s2 - u){At - 3s2) - 3t'^{ut + 2sl) + usl{5u - 2S2) + 2ts2{t^ + 2sl) 

- h'^''^^ ds{s2 - uf [ 2s^ {u{u -8t- 2s2) - t{t - 2S2)) 

+ S (3{U^ - t^) + 10S2(t^ - U^) + 8S2{4:Ut + US2 - tS2) - ut{17u + 15t) 

+ bs2{t^ - u^) + u^iu^ + 9ts2 + 8sl) - t^{f + 6ut - 23us2 + lOu^) 

- Ausl{s2 + 4:t) + 4:tsl{s2 - 2t) 

- i/f '^^ 2ut{s2 - uf - \u{3t -2u + S2) + 2sl - SS2 



3u'^{t + 5m) - 6s2(6i^ - si) + 22s2ut - -uH{17t - 47s2) 



2u^{3s2 -u)-^u (l4tsl - 3sl + tu^) | 



+ ^ I ftudsu S{Q'^ - S2- 2s) + S2(Q^ - s) - [(S2 - U){S2 -t) - SS2] 

+ fsusds \2s'^{u - 2S2) + s {u{u + bt) - 2s2{u + 2t)) 



is^iQ -s) - ut{2u - t - 3S2) + M S2 



+ fxtdsu [SS2{2S2 -U-t)+StU+ {S2 - t){2t -S2){t + U- 2^2)] 



u ^ t 



Diagrams 2{H^ + i/2)*(^3 + H^) and 2(1/5 + Hq)*{H^ + H^) for qq ^ Ym 
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The diagrams are shown in Fig. 9 of Ref. 0. The vector- vector couphng contribution for the 
projection matrix elements for the diagrams 2{Hi + H2)* {H^ + H4) and 2{H5 + Hq)*{Hj + 
Hs) are denoted by E^^yy and E^^yy, respectively, and can be obtained from the matrix 
elements D^^yy listed above: 

Eabvvis,t,u,Q'^) = E^^^yy{s,u,t,Q'^) = -D^ayy{s,u,t,Q'^) . 
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